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Abstract 

In this paper, we calculate the branching ratio and direct CP asymmetry parameter of B^r — > 
D®K^ in the framework of perturbative QCD approach based on Ut factorization. Besides 
the usual factorizable diagrams, both non-factorizable and annihilation type contributions are 
taken into account. We find that (a) the branching ratio is at the order of 10~ 5 ; (b) the tree 
annihilation diagrams and the penguin diagrams dominate the total contribution; and (c) the 
direct CP asymmetry is about 7%, which can be tested in the Large Hadron Collider beauty 
experiments (LHC-b) at CERN. 
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I. INTRODUCTION 



The rare B meson decays arouse more and more interest, since it is a good place for 
testing the Standard Model (SM), studying CP violation and looking for possible new 
physics beyond the SM. In recent years, the theoretical studies of B u d mesons have been 
done in the literature widely, which are strongly supported by the running B factories 
in KEK and Stanford Linear Accelerator Center (SLAC). B physics studies are further 
supported by the Large Hadron Collider beauty experiments (LHC-b), it is estimated 



that about 5 x 10 10 B c mesons can be produced per year at LHC 



na. 



So the studies of 



B c meson rare decays are necessary in the next a few years. 

The nonleptonic decavs of the B r mesons have been studied in previous 



S, Q, B, B, 0, H B, 0, 0, 0, 0, 3, 0, 0, 0, 0, 0, 
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employing Naive Factorization [30], QCD Factorization 3JJ, perturbative QCD approach 
(PQCD_n32j and other approaches. The theoretical status of the B c meson was reviewed 
in Ref.[l|. In this paper, we will study the B c — > DK decays in the perturbative QCD 
approach. Our theoretical formulas for the decay B c — > DK in PQCD framework are 
given in the next section. In section 1111} we give the numerical results of the branching 
ratio and CP asymmetries of B c — > DK and the form factor of B c — > D etc. At last, we 
give a short summary in section IIVI 



II. PERTURBATIVE CALCULATIONS 



For decay B c — > DK, the related effective Hamiltonian is given by 33] 

10 



H, 



cff 




E v «°K b \ c ^)° q M + c 2 (/i)o 2 V)] - v*v ts £ CiMOiQ*) , (l) 



q=u,c 



i=3 



where Ci(fi)(i = 1, ■ • • ,10) are Wilson coefficients at the renormalization scale fi and 
Oi(i = 1, ■ ■ ■ , 10) are the four quark operators 



0\ = 


{biqj)v~A{qjSi) V -A, 


oi = 


(biqi) V - A {qjSj)v-A, 


o 3 = 


(biSi) V -A J2 q (Wj)v-A, 




(biSj)v-A J2q(mdv-A, 


o 5 = 


(biSi)v-A T, q imj)v+A, 


o 6 = 


(biS^V-A T,q(mi)v+A, 


o 7 = 




8 = 


\{biSj)v-A Eg CqiW^V+A, 


o 9 = 


\(hiSi)v-A E g e q(qjqj)v-A, 


o 10 = 


-- l{b i s j ) v -. A Y. q e q(.qjqi)v-A 



(2) 
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Here i and j are SU(3) color indices; the sum over q runs over the quark fields that are 
active at the scale \i = 0(m,b), i.e., q G {u, d, s,c, b}. Operators Oi,0 2 come from tree 
level interaction, while 3 , 4 , 5 , 6 are QCD-Penguins operators and 7 , 8 , Og, O w 
come from electroweak-penguins. 

In PQCD approach, the decay amplitude can be written as: 

Amplitude - [d%d A k 2 d% Tr[C(t)$ B MQD(fo)QK(fo)H(k u k2,k 3l t)], (3) 
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where k{ are the momenta of light quarks included in each of the mesons, and Tr de- 
notes the trace over Dirac and color indices. C(t) is the Wilson coefficient results from 
the radiative corrections at short distance. $m is the wave function which describes the 
hadronization of mesons. The wave functions should be universal and channel indepen- 
dent, we can use $m which is determined by other ways. The hard part H are channel 
dependent but fortunately perturbative calculable. 

Working at the rest frame of B c meson, the momenta of the B c , D and K can be 
written in the light-cone coordinates as : 

Mb , , M B , 9 , M B , 9 

P X = -|(1,1,0 T ), P 2 = _^(r 2 ,l,0 T ), P 3 = _^(l-r 2 ,0,0 T ). (4) 

where r = M B /M Bc and we neglect the kaon's mass Mr. Denoting the light (anti-)quark 
momenta in B c , D and K as ki, k 2 and k$, respectively, we can choose: 

h = (xipf, 0, ki r ), k 2 = (0, X2P2, k 2T ), h = (x 3 p^, 0, k 3T ). (5) 

Then integration over fcf , and k% in Eq.([3]) leads to 

Amplitude ~ Jdxidx2dx 3 bidbib2db 2 b 3 db 3 

xTr[C(t)^ Bc (x 1 ,b 1 )<l> D (x2,b2)^K(x 3 ,b 3 )H(x l ,b t ,t)}e~ s ^, (6) 

where 6j is the conjugate space coordinate of k iT , and t is the largest energy scale in H. 
The exponential Sudakov factor e~ s ^ comes from the higher order radiative corrections 
to wave functions and hard amplitudes, it suppresses the soft dynamics effectively [34| and 
thus makes a perturbative calculation of the hard part H reliable. 

According to effective Hamiltonian (Op), we draw the lowest order diagrams of B c — > DK 
in Fig. [U The usual factorizable diagrams (a) and (b) give the B c — > D form factor if take 
away the Wilson coefficients. The operators Oi, 2 , 3 , 4 , 9 and O w are (V — A)(V — A) 
currents, and the sum of their contributions is given by 

F e [C] = -^=7rC F M 2 Bc / dx 2 / b l db l b2db2<t>D{x2,b 2 ) 

y£N c Jo Jo 

x {[( r 2 + 2r _ i) X2 + (2 - r)r b ]a s (t 1 a )h^(x 2 , b u b 2 ) exp[-S B {t l a ) 

- 5 (4)]C(4) + [r(2 - r)}a s (t 2 a )h { : ) (x2, h, b 2 ) exp[-S B (t 2 a ) 

- S D (t 2 a )}C(tl)}, (7) 

where r& = Mf,/M Bc , Cf = 4/3 is the group factor of the SU(3) C gauge group. The 
expressions of the meson distribution amplitudes 0a/, the Sudakov factor Sx{ti)(X = 



4 



B c , K, D), and the functions h a are given in the appendix. In above formula, the Wilson 
coefficients C(t) of the corresponding operators are process dependent. 

The operator 5 , 6 , 7 , O s have the structure of (V — A)(V + A), their amplitude is 

F e p [C] = ^^r K nC F M 2 Bc f dx 2 H ' b 1 db 1 b 2 db 2 <j> D {x 2 M) 

x {[(2r 2 - r)x 2 + (2 - r - r b - 2r 2 + Ar b r)]a s (t l a )h^\x 2 , b u b 2 ) 
x exp[-5 fl (0 - + [3r -4r 2 ]a s (t 2 )^)( a ; 2 ,6 1 ,6 2 ) 

xexp[-^(a-^(^)]C(0}, (8) 

where = M 0K /M Bc = M 2 K j\M Bc [M s + M M )]. For the non-factorizable diagrams (c) 
and (d), all three meson wave functions are involved. Using 5 function 5(bi — b 3 ), the 
integration of b\ can be preformed easily. For the (V — A)(V — A) operators the result is: 

16 f 1 f°° 

M e [C] = — 7rC F f Bc Ml c / dx 2 dx 3 / b 2 db 2 b 3 db 3 (f) D (x 2 ,b 2 )<f)£(x 3 ) 

c J J 

x {[1 -2r + (r-r 2 )x 2 - (1 - 1r 2 )x 3 \a s {t\)h^ [x 2 , x 3 , b 2 , b 3 ) 

x exp[-S- B (^) - - S^)]C(t*) + [2r - 1 + (1 - r - r 2 )x 2 

- (1 -2r 2 )x 3 ]a s (t 2 )^ 2) (a:2,a:3,&2,&3)exp[- 1 SB(t c 2 ) - S D (t 2 ) 

- S K (tl)]C(tl)}, (9) 

For the(y — A)(V + A) operators, the formula is: 

16 Z" 1 Z" 00 

M e[C] = —nC F f Bc r K M 2 Bc / ofe 2 ofe 3 / b 2 db 2 b 3 db 3 <p D (x 2 ,b 2 ) 

Jo JO 

x {[(1 + r - rx 2 - (1 + r)rr 3 )0£(a;3) + (1 - r + rx 2 - (1 + r)z 3 ) 

x ^(x 3 )]a s (^)/i«(x 2 ,x 3 ,6 2 ,6 3 ) expf-S-B^) - S^) - M*c)]C($ 

+ [(rx 2 - (1 + r):r 3 )0£(:r 3 ) + (-2r + rx 2 + (1 + r)x 3 )0^(a; 3 )] 

x a s (t 2 c )h^(x 2 ,x 3 ,b 2 ,b 3 )exp[-S B (t 2 c ) - S D (t 2 c ) - S ' K {t 2 c )]C '(t 2 )} , (10) 

Similar to (c),(d), the annihilation diagrams (e) and (f) also involve all three meson wave 
functions. Here we have two kinds of amplitudes, M a is the contribution containing the 
operator of type (V — A)(V — A), and M£ is the contribution containing the operator of 
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type (V - A)(V + A). 

16 f 1 f 00 

M 4C] = jT^C F f Bc Ml c / dx 2 dx 3 / b l db 1 b 2 db 2 (t) D {x2M) 

c J J 

x {[(1 - r - r b - r 2 x 2 + (2r 2 - l)x 3 + r 6 r 2 )</4(a; 3 ) + (2 - x 2 - x 3 - Ar b ) 
x rr K (t>K{ x z) + ( x 2 - x 3 )rr K (p r ^(x 3 ,)]a s (t 1 e )h ( ( } ) (x 2 ,x 3 ,b 1 ,b 2 )exp[-SB(t 1 e ) 
- S D (tl) - S K (tl)]C(tl) + [(r + x 2 )0^(x 3 ) + (x 2 + x 3 )rr x 0£(x 3 ) 
+ {x 2 - x 3 )rr K <p^(x 3 )]a s ( y t 2 e )h {2) (x 2 ,x 3 ,b 1 ,b 2 )exp[-S B (t 2 e ) - S D (t 2 e ) 
-S K (tl)]C(tl)}, (11) 

16 Z" 1 f°° 

M a[C] = ^C F f Bc M 2 Bc / dx 2 dx 3 / b x db x b 2 db 2 <\) D {x 2 ,b 2 ) 

Jo Jo 

x {[(-r + rx 2 - r b r)(j)K(x 3 ) + (1 + r b - r - x 3 )r K (p^(x 3 ) + (1 + r b 

- r - x 3 )r K (j) T K (x 3 )}a s (t\)h { p(x 2 ,x 3 MM) exp[-5' B (^) - S D (t\) 

- S K (t\)]C(t\) + [(r 2 - rx 2 )<f>£(x 3 ) + (x s - 2r)r K <j ) p K (x 3 ) 

+ (x 3 - 2r)r x 0^(x 3 )]a s (t 2 )^ 2) (x 2 ,x 3 ,6i,6 2 )exp[- 1 S B (te) - S D {t 2 e ) 

-S K (tl)}C(t 2 e )}, (12) 

The amplitude for the factorizable annihilation diagrams (g) and (h) result in F a (for 
(V — A){y — A) type operators) and (for (V — A)(V + A) type operators): 

pi poo 

F a [C] = \QitC F M 2 Bc / dx 2 dx 3 / b 2 db 2 b 3 db 3 (p D (x 2 ,b 2 ) 
c Jo Jo 

x {[(r 2 - l)x 2 0^(a;3) - 2rr K (l + x 2 )(j)K{x 3 )]a s (t l g )h g {x 2) x 3 ,b 2 ,b 3 ) 

x expf-S^) - S K §\]C{$) + [(-r 2 + x 3 - 2r 2 x 3 ) ( f ) £(x 3 ) 

+ (r K r + 2r K rx 3 )(pK{x 3 ) + (-r^r + 2r K rx 3 )4> T K (x 3 )} 

x a s (t 2 )/ i ,(a: 3 ,x 2 ,6 3 ,6 2 )exp[- 1 S Z) (t 2 ) - ^(t 2 )]C(t 2 )}, (13) 

pi poo 

F[[C] = 32ttC f M 2 Bc / dx 2 dx 3 / b 2 db 2 b 3 db 3 <p D (x 2l b 2 ) 
c Jo Jo 

x {[^20^(^3) + 2r K (p 1 p C (x 3 )]a s (t 1 g )h g (x 2 , x 3 ,b 2 ,b 3 ) expl-Soitl) 

- S^J)]C(^) + N40r 3 ) + r^a^^s) - r x x 3 0^(a: 3 )] 

x a s (t 2 g )h g (x 3 , x 2 , 6 3 , 62) exp[-^(t 2 ) - Sx(t 2 )]C(t 2 )}, (14) 
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From Equation (l7])- fTH|) . the total decay amplitude for 5+ — > D°K + can be written as 

A(B C - D°K + ) = f K F e [V us V: h {\Cx + C 2 ) - V t tV ts (±C 3 + C 4 + ±C 9 + C 10 )} 

- f K V t lV ts F e p [^C 5 + C 6 + l -C 7 + C 8 ] + M e [V us V* b Ct - V t tV ts (C 3 + C 9 )} 

- V*V ts M P (C 5 + C 7 ) + M^V^Cx - V*V ts (C 3 + C 9 )] 

- V^M^Cs + C 7 ) + fsMVcsV^Cx + c 2 ) - v t tv ts (~c 3 + c 4 

+ ~C 9 + C 10 )] - fB c V t * b V ts F a p [ic 5 + C 6 + ic 7 + C 8 ], (15) 

and the decay width is expressed as 

r(S+ - = ^-^(1 - r 2 )|A(5 c + - JJ°/,+)| 2 . (16) 

The decay amplitude of the charge conjugate channel B~ — > D°K~ can be obtained by 
replacing V qs V q * b to V* s V qb and V£V U to V^* in Eq.©. 

Using the unitary condition of the CKM matrix elements V* b V us + V^V^ + V^V ts = 0, 
the decay amplitude of — > D°K + in Eq. (fT5]) can be parameterized as 

a = v: b v us r u + v; b v cs T c - v t ;v ts p 

= V: b V us (T u + P)[1 + ze< 5 -% (17) 



where z = I Y,tY, a 1 1 ^ c+p I, K,/, ~ I V^,& |e n and 5 = arg( Pc+ p ). z and 5 can be calculated 

1 V* b V U s 11 T u +p 1 ' ov Tu+P' 

from PQCD. 

Similarly, the decay amplitude for B~ — ► D°K~ can be parameterized as 

= K fe Kr s (T u +p)[i+^ + ^], (is) 

and the averaged decay width for B+(B~) — ► D°K ± reads 

r(B+(B") 23°^) = -^(l-r 2 )(|A| 2 /2 + |A| 2 /2) 

G 2 M 3 

= -fz^Q- - r 2 )\V: b V us {T u + P)| 2 [l + 2ZCOS7COS5 + z 2 ](19) 
which is a function of CKM angle 7. 
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III. NUMERICAL EVALUATION 



The following parameters have been used in our numerical calculation 24|, |35j, |36j, 137 . 



M Bc = 6.286GeV, M b = 4.8GeV, M D = 1.685GeV, M 0K = 1.6GeV, 
co D = 0.2GeV, f Bc = 489MeV, f K = 160MeV, f D = 240MeV, 

0.0009, \V; b V cs \ = 0.0398. (20) 



a D = 0.3, r Bc = 0.46 x 1Q- I2 s, \V* b V us \ 



We leave the CKM phase angle 7 as a free parameter to explore the branching ratio and 




25 50 75 100 125 150 175 
y(degree) 



FIG. 2: The averaged branching ratio of B^(B C ) — ► D°K ± decay as a function of CKM angle 
7- 



CP asymmetry parameter dependence on it. The averaged branching ratio of the decay 
B+(B~) — ► D°K ± with respect to the parameter 7 is shown in Fig. [2j Since the CKM 
angle 7 is constrained as 7 = (63j^2)° 35], we can arrive from Fig. [2J 



6.5 x 10" 5 < Br{B+(B~) -> D ^) < 6.7 x 10~ 5 , for50° < 7 < 80°. 



(21) 



Our numerical analysis show that \V* b V cs T c /V* b V us T u \ = 10 and \V t * b V ts P/V* b V us T u \ = 
7.4 which mean the tree level contributions from annihilation topology and the penguin 
contributions dominate in this decay and the branching ratio is not sensitive to 7. In 
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our calculation, the only input parameters are wave functions, which represent the non- 
perturbative contributions. In all the three meson wave functions, the main uncertainty 
come from the value of in D meson wave function(see appendix). We investigate the 
branching ratio dependence on the value of ujd in Table HI By changing the value of ud 



from 0.26W 



38 



to OAGeV, we find the branching ratio of B+(B~ 



D°K ± change 



little as shown in Table HI 





oo D = 0.2GeV 


UJ D 


= 0.3GeV 


UJ D 


= OAGeV 


B+ -» D°K+ 


6.6 


7.0 


5.6 



TABLE I: Branch ratios in the unit 10 5 using 7 = 60° for different up 



The diagrams (a) and (b) in Fig. [T]give the contribution for B c — > D transition form 
factor F Bc ^ D (q 2 = M\ ~ 0), where q = P± — P2 is the momentum transfer. The sum of 
their amplitudes have been given by Eq. (JTj), so we can use PQCD approach to compute 
this form factor. Our result is F B ^ D (0) = 0.24, if u D = 0.2GeV; and F B ^ D (0) = 0.21, 
if Ur> — 0A5GeV. We can see that this form factor is not sensitive to the D meson wave 
function. In the literature, there already exist a lot of studies on B c — > D transition form 



factor 



39 



40 



41 



42 



43 



44j, we show their results in Table HT1 From which we 
find that there are large differences in these results(including ours) and eventually this 
form factor can be extracted from semi-leptonic experiments B c — > Dlv\ in the future 
LHC experiments. 





DW[27] a 


CNP[39] 


NW[40] 


IKS [41] 


KKL[12J 6 


EFG[8] 


ZH[42] 


DSV[43j 


WSI. 11 


F B ^ D {0) 


0.255 


0.13 


0.1446 


0.69 


0.32[0.29] 


0.14 


0.35 


0.075 


0.16 



a We quote the result with lo = 0.7GeV. 

^The nonbracketed(bracketed) result is evaluated in sum rules (potential model). 

TABLE II: B c — > D transition form factor at q 2 = evaluated in the literature. 



The direct CP violation Afj r P is defined as 



a dir 
A CP 



T(B: 



T(B+ -> D°K+) + T(B- 
2z sin 7 sin 5 

1 + 2z cos 7 cos 5 + z 2 



D°K~) 



[22] 
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Using Eq. (TTT|) and (TT8"j) . we compute the parameter A^ r p . The direct CP asymmetry 
A^p has a strong dependence on the CKM angle 7, as can be seen easily from Eq. (1221) 
and Fig. [31 From this figure one can see that the direct CP asymmetry at 6% — 7.6% 
for 50° < 7 < 80°. The small direct CP asymmetry is also a result of small tree level 
contribution from emission topology. 




y(degree) 

FIG. 3: Direct CP violation parameters of B^- — > D Q K^ decay as a function of CKM angle 7. 



IV. SUMMARY 

In this work, we study the branching ratio and CP asymmetry of the decays Bf — > 
D°K ± in PQCD approach. It is found that the branching ratio of Bf — > D°K ± is at the 
order of 1CT 5 . We also predict CP asymmetries in the process, which may be measured 
in the LHC-b experiments. 
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APPENDIX A: FORMULAS FOR THE CALCULATIONS USED IN THE 
TEXT 



In the appendix we present the explicit expressions of the formulas used in section II. 
First, we give the expressions of the meson distribution amplitudes 4>m- For B c meson 



wave function, we use the function as 



29J: 



$rAx) 



4N r 



( Bc + M Bc )j 5 5{x-M c /M Bc ). 



For D meson wave function, we use the function as [45]: 

3 1 
Mx,b) = -j==f D x(l - x){l + a D (l - 2x)} exp[--(uj D b) 2 ] 

We set cld = 0.3 and up = 0.2Ge^|38j in our numerical calculation. 

lA,P,T 



We use the distribution amplitude y K 



46j: 



<f> T K (x) 



2y/2N c 

Ik 



2^2N~ C 
fx 



of the K meson from Ref. 
x(l -x)[l + 0.15* + 0.405(5t 2 - 1)], 

[1 + 0.106(3t 2 - 1) - 0.148(3 - 30t 2 + 35t 4 )/8 

t[l + 0.1581(5t 2 -3)], 



where t = 1 — 2x. whose coefficients correspond to itlqk = 1.6GeV 
Sb c , Sd, Sk used in the decay amplitudes are defined as 



S Bc {t)=s(x 1 P+,b 1 )+2 



dp 
dp 



S D (t) = s(x 2 Pz, b 2 )+2 / -^ 7 (a,(/2)), 

Ji/b 2 ( i 

S K (t) = s(x 3 P+, b 3 ) + s ((l - x 3 )P+, b 3 ) 



* dp 

l/b 3 P 



(Al) 



(A2) 



(A3) 



(A4) 
(A5) 
(A6) 



where the so called Sudakov factor s(Q, b) resulting from the resummation of double 

m\ 



logarithms is given as |47l . 

s(Q,b) 

with 



/// 

1/bP 





"67 


c F ^ + 


9 


71 




e 2 7B -l 






3 71 I 


2 



m i 9) 



10 



A(a{p)) + B{a s {p)) 



-nt H — 8 In 
3 27 ; 3 



d7b 



71 



(A7) 

(A8) 
(A9) 
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Here ^ E = 0.57722 • • • is the Euler constant, rif is the active quark flavor number. 

The functions hi(i = a, c.e.g) come from the Fourier transformation of propagators of 
virtual quark and gluon in the hard part calculations. They are given as 

h^(x 2 ,h,b 2 ) = S t (x 2 )K (M B ^(l - x 2 ){r - r 2 )h) 



x [0(6 2 - h)6(r 2 b - (1 - r 2 )x 2 )I (M B Jr 2 - (1 - r 2 )x 2 h) 



x K (M B y/r 2 b -(l-r2)x 2 b 2 ) + (b, <-> b 2 )}, (A10) 
h^(x 2 ,h,b 2 ) = S t (r)K (M By /(l-x 2 )(r-r*)b 2 ) 

x [9(b 2 - b^IoiMsV^^b^KoiMsVT^h) + (6i <-> & 2 )] , (All) 

h ( c J \x 2 ,x 3 ,b 2 ,b 3 ) = 
6(b 2 - & 3 )/ (AW(l - rr 2 )(r - r 2 )& 3 )F (M B ^(1 - rr 2 )(r - r 2 )& 2 ) 

1 / K (M B F m & 3 ), for F 2 > ()\ 

+ 4:<o)' (A12) 

where Hq 1 ^) = Jq(z) + iY (z), and i^'s are defined by 

ifi) = (1 - r - x 3 + r 2 x 3 )(x 2 - 1), if 2) = (1 - x 2 )(r - r 2 - x 3 + r 2 x 3 ); (A13) 

^ } (x 2 ,x 3 ,6i,6 2 ) = 
6(b 2 - 6 1 )^Hj 1) (M i?v /(l - r 2 )x 2 x 3 6 2 )J (M BV /(l - r*)x 2 xM 

, h h A ( K (M B F e(j -)6i), for > \ 

^^^^fH^M^^for F 2 ,< J' (M4) 

where F e 0')' s are defined by 

^e(i) = - (1 - z 2 )(l - r - x 3 + r 2 x 3 ), F 2 {2) = r 2 - x 2 (x 3 - r - r 2 x 3 ); (A15) 

h g (x 2 ,x 3 ,b 2 ,b 3 ) = St(x 2 )Y H Q V ' \M B ^ (1 - r 2 )x 2 x 3 b 3 ) 

x [6(b 3 - b 2 )J {M By /{l - r*)x 2 b 2 )^HW(M B y/(l - r*)x 2 b z ) + (b 2 <- & 3 )] . (A16) 

12 



We adopt the parametrization for S t {x) contributing to the factorizable diagrams 49] , 
St(x)= r } n , , \ J K1-^)] C , c=0.3. (A17) 

y/TVL (1 + Cj 

The hard scale t^s in Eq.(|71)-(|T4l) are chosen as 

t l a = max(M BV /(l - x 2 )(r - r 2 ), M By J\r% - (1 - r 2 );r 2 |, I/61, l/6 2 ), 

t£ = max(M B v / r-r 2 ,l/6 1 ,l/6 2 ), 

£ = max(M B y^y , M B y/(r - r 2 )(l - x 2 ), l/b 2 , l/6 3 ), 

t c 2 = max(M B ^F^, M BV /(r - r 2 )(l - x 2 ), l/6 2 , l/6 3 ), 

t* = max(M B ^|F 2 1) |,M i?v /(l-r 2 )x 2 X3,l/6i,l/6 2 ), 

^ = max(M B ^|F 2 (2) |, M BV /(1 - r 2 )x 2 x 3 , l/&i, l/6 2 ), 

*J = max(M BV /(l-r 2 )x 2 , l/6 2 ,l/6 3 ), 

t 2 , = max(M BV /(l-r 2 )x 3 ,l/6 2 ,l/&3)- (A18) 

They are given as the maximum energy scale appearing in each diagram to kill the large 
logarithmic radiative corrections. 
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